(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date 
24 April 2003 (24.04.2003) 




PCT 



(10) International Publication Number 

WO 03/034780 A2 



(51) International Patent Classification 7 : H04R 3/00 

(21) International Application Number: PCT/GB 02/04605 

(22) International Filing Date: 10 October 2002 (10.10.2002) 

(25) Filing Language: English 

(26) Publication Language: English 



(30) Priority Data: 

0124352.6 



11 October 2001 (11.10.2001) GB 



(71) Applicant (for all designated States except US): 1...LIM- 
ITED [GB/GB]; St John's Innovation Centre, Cowley 
Road, Cambridge, CB4 0WS (GB). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): GOUDIE, Angus, 
Gavin [GB/GB]; 2 Canterbury Close, Cambridge CB4 
3QQ (GB). THROUGHTON, Paul, Thomas [GB/GB]; 
71 Gwydir Street, Cambridge CB1 2LG (GB). HOOLEY, 
Anthony [GB/GB]; 79 De Freville Avenue, Cambridge 
CB4 1HP (GB). 



(74) Agents: MERRYWEATHER, Colin, Henry et al.; J.A. 
Kemp & Co., 14 South Square, Gray's Inn, London WC1R 
5JJ (GB). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 

AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 
CZ, DE, DK, DM, DZ, EC, EE, ES, FT, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ, NO, NZ, OM, PH, PL, PT, RO, RU, SD, SE, SG, 
SI, SK, SL, TJ, TM, TN, TR, TT, TZ, UA, UG, US, UZ, 
VC, VN, YU, ZA, ZM, ZW. 



(84) 



Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, IE, IT, LU, MC, NL, PT, SE, SK, 
TR), OAPI patent (BF, BJ, CF, CG, CI, CM, GA, GN, GQ, 
GW, ML, MR, NE, SN, TD, TG). 



Published: 

without international search report and to be republished 
upon receipt of that report 



[Continued on next page] 



(54) Title: SIGNAL PROCESSING DEVICE FOR ACOUSTIC TRANSDUCER ARRAY 



811-n 



< 



QC 

rn 
O 



D 



o - O 
n o ° o op? o 



(J 



TJ 



o 



o o o " ° 




O a 

b° 



o 



00»Vo o "o00 

o °o 0 o 



o 



o 



o 



Q °Q 



o 



o 



c 



o 



811-1 

(57) Abstract: The invention provides transducer arrays which are capable of outputting sound beams having a relatively constant 
width, and with minimal sidelobes, across a range of frequencies. This is achieved by utilising one or more digital signal modifiers 
within the signal path between the input sound signal and the array of transducers. Variable window functions are also disclosed. 
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SIGNAL PROCESSING DEVICE FOR ACOUSTIC TRANSDUCER ARRAY 

FIELD OF THE INVENTION 
This invention relates to steerable antennae and arrays of transducers, and concerns 
5 in particular arrays of electro-acoustic transducers. 

BACKGROUND OF THE INVENTION 
Steerable or phased array antennae are well known in the art in both the 
electromagnetic and the ultrasonic acoustic fields. They are less well known in the 
10 sonic (audible) acoustic area. 

The commonly-owned published International Patent application No WO 01/23104 
describes sonic steerable or phased array antennae and their use to achieve a variety 
of effects. The application describes a method and apparatus for taking an input 
15 signal, replicating it a number of times and modifying each of the replicas before 
routing them to respective output transducers such that a desired sound field is 
created. This sound field may comprise a directed beam, focussed beam or a 
simulated origin. 

20 Control of direction and beamwidth, i.e. the steerability, of a beam is required to 
generate and steer broadband acoustic signals, such as multi-channel audio signals. 
These parameters depend on the frequency or range of frequencies of the emitted 
signal. In addition they depend on the spatial arrangement of the emitting sources. 
The spatial arrangement in turn is subject to technical constraints arising from the 

25 technical properties of the transducers employed and costs. Thus, the design of a 
functional and economically viable source of acoustic energy capable of projecting 
sound into predetermined directions, in short herein referred to as digital loudspeaker 
system or DLS, is a complex task. 

30 In WO 01/23104 the direction of a beam is controlled by delaying the output of each 
transducer across the array. Appropriate delays, which are frequency dependent, lead 
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to a constructive interference at a predetermined location of all the signals as emitted 
from the transducers of the array. 

On the other hand, the beamwidth - whether measured as the angular distance 
5 between two minima or by any other known definition - is in the simplest case a 

function of direction of the beam, its frequency and the emission area or width of the 
array of sources from which the beam emanates. For previously-described arrays, the 
beam becomes narrower with increasing frequency. With broadband signals, 
spanning a broad range of frequencies, potentially many octaves in case of audio 
10 signals, this makes it difficult to generate and steer a beam at the lowest frequency 
components of the signal. One way to overcome this problem is by extending the 
lateral dimensions of the array of the antennae. However, such larger array narrows 
the beam at high frequencies. This effect could be disadvantageous in practical 
applications such as, for example, the projection of sound. 

15 

It is therefore an object of the invention to improve the ability of an array of acoustic 
transducers to emit and steer beams of broadband sonic signal while minimizing 
mechanical and electronic components required for its implementation. 

20 It is another object of the invention to obtain an array of broadband transducers that 
emits broadband wave signal with sufficient directivity at low frequency and 
sufficient beamwidth at high frequencies. 

It is a further object of the invention to obtain an array of broadband transducers with 
25 improved steerabilty of sound beams having different travel paths before reaching a 
listener. 



30 



SUMMARY OF THE INVENTION 
In view of the above objects, the present invention provides a method and apparatus 
as claimed in the independent claims. 
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According to a first aspect of the invention, there is provided an array of electro- 
acoustic transducers capable of steering one or more beams of signal. The signal, 
being preferably an audio signal, consists of components at many different 
frequencies simultaneously present in the signal. By using appropriately configured 
5 digital signal modifiers, such as digital filters, that adjust the output response array 
for each of these different components a non-zero output can be limited to subarrays 
of the array. By broadening the borders of subarray with decreasing frequency of the 
signal components, a constant beamwidth can be achieved over a whole range of 
frequencies. 

10 

In a variant of this aspect of the invention the edge of the effective area is 
smoothened by spreading the reduction from full amplitude or gain to cut-off or zero 
output over a zone that includes at least one transducer operating at a gain level 
between those two values. The smoothened is intended to reduce the amount of 
15 energy emitted as sidelobes to the main beam or beams. 

A particularly convenient way of implementing the digital signal modifiers is as 
digital finite impulse response filters programmed to emulate a window function. 
The window function widens the area of non-zero emission with decreasing 
20 frequency, thus maintaining a constant beamwidth of the signal over a large 

frequency range. Many different window functions can be used within the scope of 
this aspect of the invention. 

It is a second aspect of the invention to introduce a physical arrangement of 
25 transducers that minimizes the number of transducers necessary to generate steerable 
beams of sonic signals. It was found that by varying the spacing between adjacent 
transducers gradually or step-wise towards the outer area of the array, the number of 
transducers could be significantly reduced in comparison with an array of equal 
width but regular spacing. Alternatively, the size of the transducers may be varied. 

30 

By considering the limitations on transducer spacing as imposed by the first aspect of 
the invention, arrays of minimal numbers of transducers can be designed, yet 
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satisfying the need to generate broadband beams of near-constant beamwidth. All of 
the above aspects are applicable to one- and two-dimensional flat or curved arrays of 
transducers. 

5 These and other aspects of inventions will be apparent from the following detailed 
description of non-limitative examples making reference to the following drawings. 



BRIEF DESCRIPTION OF THE DRAWINGS 

10 In the drawings: 

FIG. 1 illustrates an example of a multi-transducer source as described in the 

International patent application WO-0123104; 

15 FIG. 2 is a block diagram showing several signal processing stages prior to 

emission within a multi-transducer source; 

FIG. 3 is the block diagram of FIG.2 modified in accordance with an 

embodiment of the invention; 

20 

FIG. 4 is a side view illustrating the effect of the invention on the device of 

FIG. 1; 

FIG. 5A is a plot of gain window functions in accordance with a first example 
25 of the invention; 

FIG. 5B shows frequency responses of digital filters derived from the window 

functions of FIG. 5 A; 



30 FIG. 6A 



is a plot of gain window functions in accordance with a second 
example of the invention; 
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FIG. 6B shows frequency responses of digital filters derived from the window 

functions of FIG. 5 A; 

FIG. 7 is a plot of gain window functions with increased gain at lower 

5 frequencies; 

FIG. 8A illustrates a possible path pattern according to which transducers may 

be positioned within an array; 

10 FIG. 8B is an array layout generated in accordance with an example of the 
invention and the path pattern of FIG. 8 A; 

FIG. 9A shows a radial array layout of an array in accordance with an example 
of the invention; 

15 

FIG. 9B is the block diagram of FIG. 3 showing a variant in accordance with 

the array layout of FIG. 9 A; 

FIG. 10 shows a elliptical array layout of an array in accordance with a further 

20 example of the invention; and 

FIG. 11 is a flow chart illustrating steps of a method in accordance with the 

invention. 



25 

DETAILED DESCRIPTION 

Firstly there is described a known arrangement of transducers capable of steering a 
beam of sonic signal into one or more predetermined directions, also referred to as 
30 DLS (Digital Loudspeaker System). 
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The basic arrangement of FIG. 1 shows an array 10 comprising a plurality of 
spatially-distributed electro acoustic transducers 11-1 to 11 -n mounted on a common 
chassis 12 and arranged in an essentially two-dimensional array. The transducers 11 
are each ultimately connected to the same digital signal input. This input is modified 
5 and distributed to feed the transducers. Beamsteering is accomplished by adding 

delays or phase shifts to the signal to ensure a constructive interference of the signals 
stemming from the individual transducers at pre-determined locations 13, 14. For 
the purpose of the present example, these location are spots on the side or rear wall 
of a room giving sufficient reflection to redirect the sound back to a listener 15 in the 

10 room. Basic geometric calculations show that the delay is a function of the relative 
positions of the transducers of the array and the direction 0 of the locations 13, 14 
relative to the transducers 11-1 to 11 -n. Whilst determining the necessary delay or 
phase shifts is a complex task in itself, the present invention seeks to improve certain 
aspects that can be treated independently of the basic beamsteering process. For 

15 further details of the delay or phase shift aspects of the beamsteering, reference is 
made for example to the published International patent application WO-0123104, 
fully incorporated herein. 

Whereas the calculation of delay and phase shifts is a known mathematical problem, 
20 the electric and electronic circuitry necessary to modify the signal such as to feed 
appropriately delayed replicas of the signal to each transducer of the array can vary 
widely and is of course subject to technological advances in the field of signal 
processing. The components of FIG. 2, as referred to in greater detail below, are 
therefore considered as being highly interchangeable with other components having 
25 the same digital signal processing capabilities. 

In FIG. 2, audio source data is received by the DLS via inputs 21 as either an optical 
or coaxial digital data stream in the S/PDIF or any other known audio data format. 
The data may contain simple two-channel-stereo signal or modern compressed and 
30 encoded multi-channel sound reproductions such as Dolby Digital 1 ™ 5. 1 or DTS tm 
sound. Multi-channel inputs 21 are first decoded and decompressed using digital 
signal processing devices and firmware 22 designed to handle these proprietary 
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acoustic data formats. Their output is fed into three pairs of channels 23. In turn, the 
channel pairs provide the input to a multi-channel sample rate converter 24 for 
conversion to a standard sample rate and bit length. The outputs of the sample-rate- 
converter stage 24 are combined into a single high-speed serial signal comprising all 
5 six channels. In case of a conventional stereo input, only two of these may contain 
valid data. 

The serialized data enters Digital Signal Processing (DSP) unit 25 to further process 
the data. The unit comprises a pair of commercially available Texas Instruments 
10 TMS320C6701 DSPs running at 133MHz and performing the majority of 
calculations in floating point format. 

The first DSP performs filtering to compensate for the irregularities in the frequency 
response of the transducers used. It provides four-times over-sampling and 
15 interpolation to remove high-frequency content generated by the oversampling 
process. 

The second DSP performs quantization and noise shaping to reduce the word length 
to nine bits at a sample rate of 195kHz. 

20 

The output from the second DSP is distributed in parallel using bus 251 to eleven 
commercially available Xilinx XCV200 field programmable gate arrays (FPGAs) 26. 
The gate arrays apply a unique time delay for each channel and for each transducer. 
Their output is a number of different versions or replicas of the input, the number 

25 being equal to the number of transducers times the number of channels. As the 
number of transducers 211-1 to 211-n in this example is 132, several hundred 
different versions or replicas of the input are generated at this stage. The individual 
versions of the channels are summed at adders 27-1 to 27-n for each transducer and 
passed to pulse width modulators (PWM) 28-1 to 28-n. Each pulse width modulator 

30 drives a class-D output stage 29-1 to 29-n whose supply voltage can be adjusted to 
control the output power to the transducers 211-1 to 211-n. 



WO 03/034780 



PCT/GB02/04605 



8 

System initialisation is under the control of a micro-controller 291. Once initialised 
the micro-controller is used to take direction and volume adjustment commands from 
the user via an infrared remote controller (not shown), display them on the system 
display, and pass them to the third DSP 292. 

5 

The third DSP in the system is used to calculate the required time delay for each 
channel on each transducer to be able to steer, for example, each channel into a 
different direction. For example, a first pair of channels can be directed to the right 
and left side-walls (relative to the position of the DLS) of a room while a second pair 
10 is directed to the right and left of the rear- wall to generate a surround sound. The 
delay requirements, thus established, are distributed to the FPGAs 26 over the same 
parallel bus 251 as the data samples. Most of the above steps are described in more 
detail in WO-0123104. 



15 Referring now to a first embodiment of the invention as shown in FIG. 3, an 

additional filtering process 31 is added to the signal path of FIG. 2. In should be 
noted that in order put emphasis on the changes introduced by the present invention 
the same reference numerals and characters designate like parts in FIGs. 2 and 3, 
respectively. 

20 

In FIG. 3, digital filters 31-1 to 31-n are applied after the signals have been separated 
according to channel and added. The output of the digital filter stage is sent to the 
PCM stage 28-1 to 28-n of each of the transducers 211-1 to 211-n. The digital filters 
31-1 to 31-n can be implemented by separate DSPs or gate arrays, or, in fact, may 
25 just be included into other signal processing devices 25, 26. 

As the physical implementation of the digital filters may vary in accordance with the 
electronic components used to build the DLS, the filters are better described in terms 
of their desired response or effect on the signal. 

30 

The filters are designed to control or modify the output of the transducers depending 
on the frequency of the signal to be emitted. Within a frequency range of 500 Hz to 
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10 kHz the filters 31-1 to 31 -n seek to maintain an approximately constant 
beamwidth. This is done in practical terms by imposing frequency dependent 
windows onto the output amplitude of the transducers 211-1 to 211-n of the array. 
Hence, the new filters reduce the gain of transducers depending on their relative 
5 position within the array and on the frequency content of the signal to be emitted. 

In the following section, making reference to FIGs. 4 to 6, this embodiment of the 
invention and further variants thereof will be described in more detail. 

10 In FIG. 4, there is illustrated the effect a device in accordance with an embodiment of 
the invention has on the operation of an array 10 of transducers 11-1 to 11-n. Again, 
the numerals used in FIG. 4 are equal to those used in FIG. 1 for equal or equivalent 
elements. 

15 The two-dimensional plots 41, 42, 43 shown in FIG. 4 illustrate the output gain 

applied to the transducers of the array at three different frequencies fl, f2 and f3 in 
order of increasing frequency. The transducer array defines a plane having a point of 
origin 441 or zero point located at the centre of the array 10. Perpendicular to the 
plane as defined by the array, there is shown a virtual axis 44 representing the gain of 

20 the emitted signals. An arbitrary albeit high attenuation is defined as the cut off level 
and drawn to coincide with the plane of the transducer array. Thus, the curves 411, 
421, 431, representing the cut-off level for signal content having a frequency fl, f2 
and f3, respectively, indicate which of the transducers of the array 10 contribute to 
the emission: Transducers positioned within the boundary set by curve 411 

25 contribute to the emission of signal having the frequency fl, transducers positioned 
within the boundary set by curve 421 contribute to the emission of signal having the 
frequency £2, and so forth. Transducers located outside the respective boundaries are 
operated at cut-off gain or below. The area enclosed by curves 411, 421, 431 are 
three representatives of what in the following is referred to as the effective emission 

30 area of the array at a given frequency f. 
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10 

It is now a purpose of the invention to control the effective emission area within 
limits mainly set by frequency and physical dimensions of the array as a means to set 
or select a frequency independent beamwidth. By varying the effective area as a 
function of the frequency this selected beamwidth can be held at constant or near 
constant value over a broad range of frequencies, typically an octave or more. To 
this end, use is made of the functional relation between beamwidth and the linear 
dimensions of the effective emission area. In the simplest case of a one-dimensional 
array of (infinitely small) sources this functional relation can be represented by 
formula [1]: 



2f sin 0 BW 



wherein l eff is the effective half length of the array at the frequency f for a given 
15 beamwidth 0 B w (given as the angle between the two minima limiting the main 
beam). The constant c is the speed of sound in air. 

Thus, by selecting a beamwidth 0 B w adapted to the specific environment in which the 
invention is sought to be implemented, the signal processing devices 31-1 to 31-n of 
20 FIG. 3 can be programmed to reduce the output of the transducer in a frequency- 
dependent manner to generate an effective emission area in accordance with formula 
[1]. 

However, the application of [1] assumes a sudden drop of the emitted signal from 
25 full to zero signal amplitude at the edge of the effective area. In the context of FIG.4, 
the attenuation plots 41, 42, 43 would depict, instead of a smooth increase to full 
signal strength, a single step to full strength at the boundary curves 411, 421, 431, 
equivalent to the application of a rectangular window. However, introducing a sharp 
edge into the emission area is likely to cause an undesirably high amount of energy 
30 to be emitted in side-lobes, i.e., less directed sound. Therefore, there a more preferred 
variants of the invention to be described in the following, which variants spread the 
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edge zone over a broader transition zone surrounding the effective emission area. 
Within this area the transducers are controlled such that their gain is gradually 
reduced to zero depending on their radial distance from a centre of the array. In FIG. 
4, the transition zone is illustrated in a disproportional manner leading to very 
5 pointed attenuation profiles or windows. In practice any known window function 
with tapering edges can be applied to create an effective emission area with a 
transition zone at the edge. 



The choice of the window function is determined by a compromise between desired 
10 beam width and sidelobe level. Suitable window functions include the Hann window, 
which can be represented by formula [2-1] 



2 7ir I I 

[2-1] w(r) = COS ( ),if r < a (and zero otherwise). 



2a 



15 For the Hann window having a relation linking the effective half length l e ff of the 
window with frequency at a given beam width 0bw is: 



[2-2]l eff 



f sin 0 BW 



20 Another applicable window is the cos window represented by 



7ir | | 

[3-1] W(r) = COS ( ), if r < a (and zero otherwise). 

2a 



For the cos window, the equivalent of relation [2-2] can be written as 

25 

3c 



[3-2] 1 



eff 



4f sin 9 BW 
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Other applicable window functions include Hamming-, Kaiser- or Chebyshev-type 
windows or windows of the sin(x)/x type (which become Bessel functions in two 
dimensions), all of which are widely documented. 

5 

Application of such window functions leads to a modified relation [1], [2-2] and [3- 
2] between frequency and effective array length. 

The use of these tapered window functions broadens the effective length l e ff 
10 compared to formula [1] which represents a box-car window. However, the general 
characteristic of [1] holds, i.e. to maintain a constant beamwidth the effective 
emission area needs to be decreased with increasing frequency and vice versa. 

After selection of a suitable window function, a set of desired filter responses can be 
15 derived from it, as shown when referring to FIGs. 5 A and 5B below. Using standard 
design tools the desired filter response can then be converted into filter coefficients 
that implement the filter in the digital domain. A known method to derive from the 
filter response the filter coefficients is for example using an inverse Fourier 
transform. Known mathematical or engineering programs, such as MATLAB™, are 
20 readily capable of performing the necessary conversion steps. The filters of this 
embodiment are linear phase finite impulse response filter, as it is regarded as 
beneficial to maintain phase relationships and delays introduced through the beam 
steering process. 

25 Alternative filter architectures, such as infinite impulse response filters with all-pass 
phase correction stages can be used. 

Independent of the filter architecture, it is possible to perform the complete signal 
processing, including the control process of the present invention and known 
30 beamsteering methods within a single digital signal processing step. 
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Again, many of the filter parameters (e.g. length of the filter, gain etc) are subject to 
constraints determined by the available electrical and electronic components. For an 
audio system the constraints are further determined by the necessity to shape the 
signal in real-time at audio frequencies, i.e. between 20Hz and 20 kHz. 

5 

As stated before, the effective emission area decreases with increasing frequencies, 
leaving fewer and fewer transducers to contribute to the output signal. Conversely, as 
the frequency decreases, the area increases. This general property leads to further 
advantageous modification of the window shape and thus the filter design. 

10 

Firstly, as the width of the window shrinks towards higher frequencies and taking 
further into account the finite width of any transducer, eventually only a transducer 
placed at the very centre of the array reproduces the highest frequencies. These 
frequencies are, therefore, not steered at all. 

15 

By setting a minimum window width, it can be ensured that a sufficient number of 
transducers are within the window radius at the cut-off level to give the signal some 
steerability. Applying a minimum window width causes the beam to further narrow 
at higher frequencies, but, depending on the application, that may be preferable to 
20 having no directivity at all. 

At the low frequency limit, i.e., as the window reaches the physical width of the 
array, several different window designs can be applied. Each of the designs has 
advantages and weaknesses with respect to different aspects of the sound emission 
25 process. 

In the example of the present invention as illustrated by FIG. 5A, a minimum and a 
maximum window are set to accommodate for the physical limits of the array. The 
plots of FIG. 5 A are one-dimensional graphs of a Hamming-type window function 
30 showing amplification or gain (in dB) factor versus radial distance (in meters) from 
the centre. The window function is plotted at ten different frequency values ranging 
from 10 kHz to 40 Hz. However, due to the implementation of a minimum and 
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maximum window, the plots for 10 and 20 kHz at the high frequency end and for 
600, 300, 150, 80 and 40 Hz at the high frequency end are identical. The plots for 5 
kHz and 2.5 kHz and 1.2 KHz are shown as separate curves. The cut-off is set at an 
attenuation of -22 dB, the lower bound of the Hamming window. The limiting 
5 curves at 10 KHz and 600 Hz, respectively, represent the high and low frequency end 
to ensure a minimum width and a maximum width of the window. In the example, 
curve 10 Khz applies to all frequencies above 10 kHz, thus ensuring that steerability 
is maintained above this frequency. Curve 600Hz applies to all frequencies below 
600 Hz avoiding a sudden change in low frequency signal level at the edge of the 
10 array. This variant suppresses sidelobes, but at the expense of a low utilisation of the 
transducers at the fringe of the array. 

Having determined the desired shape of the windows, digital filters can be derived 
therefrom. 

15 

To derive a digital filter for transducer located for example at position R = 0.64 m, a 
frequency response characterizing the filter is obtained (conceptionally) by 
registering the attenuation values against the frequency values taking vertical section 
at position R through the window function of FIG. 5 A. As can be seen, the cut-off 
20 frequency at R=0.64m is below 2.5 kHz. Towards lower frequencies the filter gain 
increases rapidly until curve for 600 Hz is reached. The corresponding attenuation 
value of -ldB is maintained by the filter for all frequencies below 600 Hz. 

In FIG. 5B, there are shown filter frequency responses for transducer positions of 
25 1.28 m, 0.64 m as described above, 0.32 m, 0.16 m, 0.08 m, 0.04 m, 0.02m and 0.01 
m, respectively. The distances are measured as radial distance from the centre of the 
array. 

It should be noted that the use of discretely spaced transducers implies that the above 
30 continuous treatment of the window function is only a rough approximation. 

However the effects of the discrete nature of the transducers are equivalent to those 
arising from the approximation of an integral by a Riemann sum and can be equally 
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compensated for. For example, when calculating the filter response from a given 
window function, the discrete spacing of the transducer can be accommodated for by 
the trapezoid rule. Application of the trapezoid rule weights the window function at 
any discrete point with a factor proportional to the distance between adjacent 
5 transducer positions. Higher order approximations, such as polynomial based or 
other, can also be used. 

Given a numerical representation of the window functions or an equivalent frequency 
response of a digital filter and applying it to the above mentioned filter design tool 
10 derives filter coefficients that can be loaded into the digital filters shown in FIG. 3. 
The filter coefficients derived by the above steps vary continuously over the range of 
frequency and radial locations that are important to the application in questions. 

In FIG. 5, a limiting curve at 600 Hz has been introduced to apply to all frequencies 
15 below the frequency at which the window width and thus the effective emission are 
would exceed the limits of the physical array. Effectively, this imposes a tapered or 
smooth emission at the edge of the array for the full frequency range or bandwidth of 
the signal. However, other implementations are possible that increase the usage made 
of the outer transducers of the array. 

20 

In the example illustrated by FIGs. 6A and 6B, the effective emission array is 
allowed to grow beyond the physical limits of the array. In FIG. 6A, a number of the 
one-dimensional graphs of the window function show amplification or gain (in dB) 
factor versus radial distance (in meters) from the centre for 10 kHz, 5 kHz, 2.5 kHz, 

25 1 .2 kHz, 600 Hz, 300 Hz, 150 Hz, 80 Hz and 40 Hz, respectively. A minimum 

window is imposed. However, the window functions of FIG. 6 A have a finite output 
level beyond 2 metres, whereas the all windows of FIG. 5 A drop to zero at this 
radius or even smaller radial positions. In terms of output of the transducers, a 
comparison of FIGs. 5B and 6B, both showing the response function at the same set 

30 of radial positions, demonstrates a generally higher output level in the response 

functions of FIG. 6B at low frequencies. However, the general level of the output is 
increased at the cost of introducing a step change in output level at the edge of the 
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array. This step increases with decreasing frequency and, in turn, may result in more 
low-frequency energy being emitted into sidelobes. 



Another approach to address the finite length of the array is to use a family of 
5 window functions: As the frequency of the first window function reaches a value at 
which the function essentially covers the whole width of the array, i.e. each 
transducer is being used, windows of the same width but with increasing average 
value could be used to improve the low frequency power output without introducing 
discontinuities. In the example as illustrated by FIG. 7, a cos x window function is 
10 used, wherein the power x equals 2 for all frequencies where the window is equal to 
or smaller than the array width. As the window reaches the limits of the array and the 
frequency is decreased further, ever-smaller values of x are selected for the window 
function. As shown in FIG. 7, this increases the amplitude or gain levels while the 
maintaining the width of the window. 

15 

According to the above embodiments, each transducer has a separate filter depending 
on its radial position. However, it is possible to exploit rotational symmetry or 
approximate rotational symmetry to reduce the number of filters. In cases where a 
number of transducers share a radial position having different angular coordinates, 
20 e.g., are arranged on a circle, these transducers will require the same low-pass 

filtering, so their input signals can advantageously be multiplexed through common 
filters. 



Also, different beamwidths can be applied to different channels of the digital 
25 loudspeaker system. Audio channels projected at more distant walls may require a 

minimal beamwidth whereas channels projected at surfaces closer to the DLS may be 
advantageously operated employing a broader beamwidth. By chosing different 
beamwidth 8 B w in the formulae [1], [2-2], [3-2] or any equivalent relation, different 
sets of windows and, hence, different sets of filters are generated, which in turn can 
30 be applied to these different channels. 
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It will be appreciated by a skilled person from the above description that the gist of 
the above described embodiments of the invention is to give the user a high degree of 
control of the output characteristic of the DLS. While being applicable to any array 
of transducers, in particular the known regularly spaced array of transducers as 
5 shown in FIG. 1, the invention seeks to take advantage of the improved control by 
introducing arrays with irregular spacing between the transducers. From the 
description below, it will be appreciated that the irregular array designs as proposed 
by the present invention share a less density of transducers at the outer fringes of the 
array. In other words, the spacing between the transducers increases with distance 
10 from the centre of the array. An extremely important advantage of this aspect of the 
present invention is to significantly reducing the number of transducers required to 
generate a steerable broadband signal beam compared to known array designs. 



To prevent sidelobes caused by spatial aliasing, the maximum spacing between array 
15 elements must be less than some fraction of the wavelength of the highest frequency 
of interest that they are emitting. This fraction is best chosen to be in the range of 
0.25 to 0.5. For broadband arrays, whose size is determined by the lowest frequency 
of interest, this constraint, when combined with a uniform spacing can result in a 
very large number of transducers. However, the maximum allowable spacing, is 
20 proportional to the highest frequency being reproduced at any point within the array. 
Since with the above window design only the central array elements reproduce the 
highest frequencies, this is the only area that needs the highest transducer density, 
and elements can become gradually wider spaced towards the edges of the array. 



25 In a further variant of the array layout, larger transducers are advantageously used 
where the spacing of individual transducers becomes wider, i.e. towards the outside 
of the array. Larger transducers are more efficient at producing low sound 
frequencies. However, ready usage of large transducers is restricted by a technical 
phenomenon generally referred to as "high-frequency beaming". High-frequency 

30 beaming is the (undesired) directional radiation from a pistonic transducer arising 

when the diameter of the transducer is of the order of the wavelength or larger. In the 
present example, however, any transducer which is small enough to satisfy the 
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maximum allowable spacing is also small enough to have negligible beaming effects, 
as its diameter is much less than a wavelength. 

For broadband arrays, it may be advantageous to use two, three or more sizes of 
5 transducer. Where several dissimilar types of transducer are used together in an 
array, it may be necessary to use filters to compensate for their differing phase 
responses. 

Although ideally the whole array is used to reproduce the lowest frequencies, a small 
10 area at the centre of the array (i.e. the small and densely packed transducers) can be 
excluded by appropriate band filtering, e.g., by placing a high-pass filter in the signal 
path transmitting the signal to these central transducers. Or, the frequency response, 
more specifically a poor low-frequency response of the transducer can be directly 
exploited to achieve a similar effect. The steerability of the beam is largely not 
15 adversely affected by such barring of low-frequency output from the central 
transducers, if the central area has a diameter that is a fraction of the signal 
wavelength in question. This idea can be generalised to encompass several types of 
transducers, each with a different low-frequency cut-off. 

20 Since the filters for the densely packed array transducers in the central area of the 
array have high cut-off frequencies and a smooth response at low frequencies, 
relatively short finite-impulse-response (FIR) filters can be used. For transducers 
closer to the fringe of the array, the cut-off frequencies are much lower, so usually 
longer filters are used. In the above embodiment, however, these outer transducers do 

25 not emit the high frequency content of the signal. Therefore, it is readily feasible to 
use multirate signal processing and downsample the signal to be emitted by the outer 
transducers to a fraction of the original sample rate, allowing the use of shorter filters 
while maintaining the degree of control. 

30 In variants using a non-uniform distribution of transducers within the array, it is 

further found to be advantageous to ensure a uniform output per unit area of the array 
prior to the application of windowed emission. This is conveniently done be scaling 
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the output of each transducer by an appropriate factor. This factor is for example 
inversely proportional to the output per unit area at the location of the transducer. 
Having a uniform power output facilitates the application of the above aspects of the 
present invention. However as above, the general nature of digital signal processing 
5 allows folding this scaling process into the general filtering process resulting in one 
set of filters. 

There are many ways to design arrays that conform to the above constraints. The best 
approach may be to use a numerical optimisation technique. However, in the 
10 following section a deterministic but sub-optimal approach is described that has the 
advantage of producing visually pleasing layouts. 

According to this example, a grid is formed covering the dimensions of the proposed 
array. Although a uniform grid could be used, since placement accuracy becomes 
15 less important with lower frequency transducers, an irregular spacing with high 
density in the middle of the array is more efficient. 

The following parameters are given at the onset of the design process: 

20 X 7 Y Dimensions of the array 

m Minimum practicable spacing for the transducers 



25 



f_max 



Alpha 



Beta 



c 



(one type only for simplicity) 

Maximum acceptable fraction of a wavelength transducer spacing 

The desired ratio of array width to wavelength 

The maximum frequency to be reproduced by the array 

Speed of sound 



Following a square spiral path over the grid, starting in the centre, expanding to 
cover the whole array, at each location: 



30 



o Evaluate the distance r of the current location from the centre 



o Evaluate the cut-off frequency/^c = mm((Beta*c)l(2*r) 7 f_max) 
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o Evaluate the minimum permissible transducer spacing 

s = c* Alpha/ f_c 
o Evaluate the practicable spacing s_p = max(5,m) 
o Evaluate the distance to the centre of the nearest already-placed 

transducer, s_m 
o if s jn > s _p, place a transducer here 



Beta can have different values horizontally and vertically, to allow for elliptical 
10 beams. For DLS projectors, this cam be used to improve for example the horizontal 
steerability for a given number of array elements or transducers. 

To ensure the greatest low-frequency directivity for a given array size, transducers 
can be manually placed at the extremities of the array when initialising the above 
15 algorithm. When then executing the algorithm, the position of the other transducers 
is calculated taking any initially placed transducers into account. 

Grid locations on the array need not to be visited in a spiral sequence. Following 
other paths results in arrays with different properties. Good symmetry, resulting in a 
20 visually appealing product, can be achieved by following a path as shown (for a very 
small grid) in FIG. 8A where the grid points are visited in the sequence of the 
numerals assigned to it. 

FIG. 8B shows an array designed using this method, with a greater value for Beta 
25 horizontally than vertically. Transducers 811-1 to 811-u are placed such that the 
above described constraints are met. Also, the transducer vary in size, with smaller 
diameter transducers positioned at the centre of the array. 

An alternative approach to designing layouts of the transducer array is to use 
30 concentric rings of transducers. Starting with one transducer in the middle of the 

array, rings are added with the increase in ring radius and the number of elements in 
the ring chosen to satisfy the maximum permissible transducer spacing, as evaluated 
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in the previous array layout algorithm. FIG. 9A shows an array generated by this 
method with transducers arranged in six concentric rings 911-2 to 911-7 with one 
transducer 911-1 located at the centre. Transducers at the two outer rings 911-6, 911- 
7 are of larger diameter than those in the centre. 

5 

FIG. 9B is a block diagram of a possible implementation of the signal processing 
required for such an ordered array. An audio signal input 921 enters high-pass filter 
922 that removes low frequency components of the signal from the part of the signal 
to be emitted by the smaller central transducers. A stage 923 removes high frequency 

10 content from the part of the signal to be emitted by the larger transducers 911-6, 911- 
7 at the outer fringes of the array and resamples the remaining signal at a lower 
sample rate. It should be noted that this and later resampling does not cause a loss or 
deterioration of the signal as the later filtering stages that implement the effective 
emission area ensure that the outer transducers do not contribute to the high- 

15 frequency components of the signal. 

Signal correction filters, 93-2 compensate for the differing amplitude and phase 
responses of the smaller and larger transducers. 

20 As the single centre transducer 911-1 will always emit all high-frequency 

components, the signal of the compensation stage 93-1 enters directly into a digital 
signal processing and delay adding stage 96-1 that is equivalent to a combination of 
stages 26, 27, 28 and 29 of FIG. 2. This stage provides the appropriate delays, 
modulation etc. necessary to control and drive the transducer for a beam steering 

25 operation of the DLS. In the signal path to the innermost ring of small transducers, 
there is a first filter 931-1 implementing a window function in accordance with the 
invention. In the signal path to the wider ring of small transducers, the signal passes 
through a further downsampling stage 924 before entering into a second filter 931-2 
to implement the window function. Similar stages of filtering 931-3 to 931-5 and 

30 downsampling 925 towards transducers located further away from the centre are 
present in the signal path to the large transducers. 
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In accordance with the variant, the each of the filters 931-1 to 931-5 are shared 
between all the transducers within one ring. And, thus, the number of computational 
operations on the signals is significantly reduced by effectively exploiting the 
symmetry of the layout. This contrasts with the scattered arrays described in FIG. 
5 8B, which may have only 2 or 4 transducers sharing the same filter. 

It is possible to extend the ordered array approach to use non-circular 'rings'. This 
corresponds to the use of non-circular window functions. Using differing Beta values 
on each axis (as in FIG 8B) corresponds to an elliptical window function. 

10 

This can be implemented in an ordered array by using elliptical rings, as illustrated in 
FIG. 10. Placing transducers around an ellipse with equal chord distances is non- 
trivial mathematically, but can be accomplished numerically using known 
algorithms, such as the binary chop algorithm. 

15 

In the example illustrated by FIG. 10, transducers 1 1 1-1 to 1 1 1-n are shown. The 
horizontal Beta as referred to above is greater than the vertical one. The maximum 
permissible transducer spacing limit is just met around each ellipse and between the 
ellipses on the horizontal axis. However, the spacing between the ellipses is closer 
20 than necessary to meet this limit at all other angles. Hence, the design uses more 
transducers than would be necessary using a non-ordered layout with the same 
parameters. It may, nevertheless, be the preferred solution, due to reduced DSP 
requirements. This approach can be further generalised to other shaped 'rings', such 
as rectangles and hexagons with correspondingly shape windows. 

25 

In FIG 11, three steps 112, 113 and 114 are shown that illustrate the sequence of 
operational steps in accordance with an example of the invention. After choosing a 
desired beamwidth or a plurality of beamwidths, a window function is selected to 
control the emission characteristics, i.e, the effective emission area in accordance to 
30 the formulae [1], [2-2], [3-2] or other similar functions. Then, filters are designed and 
programmed to impose the window function onto the outputs of the transducers of 
the array. In operation the filters ensure that the emission is correctly widened or 
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narrowed to ensure a constant beamwidths or constant beamwidths over the range of 
frequencies present in the signal to be emitted. 

The above steps can be applied to transducers arrays of any layout. However, the 
5 layout may be optimized in accordance with further steps described hereinabove. 



The above-described methods for designing an array layout based on a window 
function produce an array that, when used with the corresponding filters, just meet 
the required condition for Alpha across the range of frequencies, thus avoiding 

10 spatial aliasing. When using smaller windows that decrease the effective emission 
area below its optimal size, beams with wider beamwidth are generated. As stated 
above, this effect when properly incorporated into the digital signal processing 
architecture can be used to control the beamwidth on a channel-by-channel basis. 
Thus, the window function used for the array layout determines a lower limit for the 

15 beamwidth, as attempting to generate a narrower beam will lead to spatial aliasing. 

The above refers to a beam at a given direction, more specifically to a direction 
perpendicular to the array. This is the direction of minimum beamwidth for a given 
array and the beams in other directions are broader. However, the methods presented 
20 above can also be used to maintain a constant beamwidth for beams in different 
directions by reducing the effective emission areas the perpendicular direction, the 
beamwidth can be held constant at a value that is sub-optimal in perpendicular 
direction but offers a constant value over most of the desired directions. 



25 
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CLAIMS 

1 . An array of transducers comprising 

5 a plurality of electro-acoustic transducers being positioned within an outer 

array boundary; 

a digital signal path between an input and said transducers for broadband 

signals with signal components in a range of frequencies; and 

one or more digital signal modifiers located within the signal path between 

10 said input and said transducers and capable of controlling outputs of said 

transducers, said one or more digital signal modifiers being adapted to 
confine outputs generated in response to said signal components to a subset 
of said transducers positioned within a subarray of said array having an 
outer subarray boundary lying within said outer array boundary, wherein 

15 said outer subarray boundary is widened quasi-continuously with 

decreasing frequency of said signal components. 

2. The array of claim 1 wherein the one or more digital signal modifiers are 
adapted to gradually reduce the output of transducers positioned within a 

20 transitional zone of the subarray from full output to effectively zero output. 

3. The array of claim 1 or 2, wherein the one or more digital signal modifiers 
are adapted to reduce the output of at least one transducer positioned within 
the transitional zone of the subarray to an amplitude level having a value 

25 below full amplitude level and effectively zero amplitude level. 

4. The array of any one of claims 1 to 3, wherein the one or more digital 
signal modifiers are adapted to widen the outer subarray boundary towards 
the outer array boundary to effectively maintain a beamwidth at a pre- 

30 selected and constant or near-constant value over the range of frequencies. 
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5. The array of any one of claims 1 to 4, having digital processors adapted to 
arrange the signal into two or more channels, said channels having different 
travel lengths to a given location, wherein the one or more digital signal 
modifiers are adapted to maintain a different beamwidth for each of said 

5 two or more channels. 

6. The array of any one of claims 1 to 5, wherein the digital signal modifier is 
a finite digital filter. 

10 7. The array of any one of claims 1 to 6, comprising further digital signal 

processors to steer one or more beams of said signal into predetermined 
directions. 

8. An array of transducers comprising 

15 a plurality of electro-acoustic transducers being positioned within an outer 

array boundary; 

a digital signal path between an input and said transducers for broadband 
signals with signal components in a range of frequencies; and 
one or more digital signal modifiers located within the signal path between 
20 said input and said transducers and capable of controlling outputs of said 

transducers, said one or more digital signal modifiers being adapted to 
impose a frequency dependent spatial gain window onto the array of 
transducers. 

25 9. The array of claim 8, wherein the width of the spatial gain window is a 

function of the frequency of the signal components. 

10. The array of claim 8 or 9, wherein the window function has a tapered edge 
at which the gain is gradually reduced with increasing window radius. 

30 
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The array of claim 8, 9 or 10, wherein the window function is independent 
of frequency for all frequencies above a higher threshold frequency within 
the range of frequencies. 

The array of claim 8, 9, 10 or 1 1, wherein the window function is 
independent of frequency for all frequencies below a lower threshold 
frequency within the range of frequencies. 

The array of any one of claims 8 to 12, wherein one or more different 
window functions are imposed for all frequencies below a lower threshold 
frequency within the range of frequencies. 

An array of transducers for creating a wavefield comprising 

a plurality of electro-acoustic transducers emitting acoustic wave signals 

and being positioned within an outer array boundary; and 

a digital signal path between an input and said transducers for broadband 

signals including signals within at least one range of frequencies, 

wherein the spacing between transducers is non-uniform within at least a 

subarray of said array. 

The array of claim 14, wherein the average distance between adjacent 
transducers increases with increasing distance of said transducers from a 
centre of the array. 

The array of claim 14 or 15, wherein transducers of a first size are 
positioned in a central subarray of the array and transducers of a second 
larger size are positioned outside said central subarray. 

The array of claim 14, 15 or 16, wherein a group of transducers are 
connected to the same one or more digital signal modifiers. 

An array of transducers comprising 
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a plurality of electro-acoustic transducers being positioned within an outer 
array boundary; 

a digital signal path between an input and said transducers for broadband 
signals with signal components in a range of frequencies; and 
5 one or more digital signal modifiers located within the signal path between 

said input and said transducers and capable of controlling outputs of said 
transducers, said one or more digital signal modifiers being adapted to 
confine outputs generated in response to said signal components to a subset 
of said transducers positioned within a subarray of said array having an 
10 outer subarray boundary lying within said outer array boundary, wherein 

said outer subarray boundary is widened quasi-continuously with 
decreasing frequency of said signal components and wherein the spacing 
between transducers is non-uniform within at least said subarray. 

15 19. A method of operating an array of electro-acoustic transducers comprising 
the steps of controlling the outputs of said transducer such that outputs 
generated in response to signal components having a range of frequencies 
are confined to a subset of said transducers positioned within a subarray of 
said array having an outer subarray boundary lying within said outer array 

20 boundary and widening said outer subarray boundary quasi-continuously 

with decreasing frequency of said signal components. 

20. The method of claim 19 comprising the step of using a frequency- 
dependent spatial gain window function to confine the outputs. 

25 

21. The method of claim 19 or 20, comprising the step of widening the outer 
subarray boundaries such that a constant or near-constant beamwidth over 
the range of frequencies is maintained. 

30 22. Sound system to reproduce a multi-channel surround sound signal including 
at least one rear channel, said system including an array of transducers 
comprising 



WO 03/034780 



PCT/GB02/04605 



28 

a plurality of electro-acoustic transducers being positioned within an outer 
array boundary; 

a digital signal path between an input and said transducers for broadband 
signals with signal components in a range of frequencies; and 
5 one or more digital signal modifiers located within the signal path between 

said input and said transducers and capable of controlling outputs of said 
transducers, said one or more digital signal modifiers being adapted to 
confine outputs generated in response to said signal components to a subset 
of said transducers positioned within a subarray of said array having an 
10 outer subarray boundary lying within said outer array boundary, wherein 

said outer subarray boundary is widened quasi-continuously with 
decreasing frequency of said signal components. 



23. The sound system of claim 22, wherein the one or more digital signal 
15 modifiers are adapted to widen the outer subarray boundary towards the 

outer array boundary to effectively maintain a beamwidth at a pre-selected 
and constant or near-constant value over the range of frequencies. 



24. The sound system of claim 22 or 23, having digital processors adapted to 
20 arrange the signal into two or more channels, including the at least one rear 

channel, said channels having different travel lengths to a given location, 
wherein the one or more digital signal modifiers are adapted to maintain a 
different beamwidth for each of said two or more channels. 



25 25. The sound system of claim 22, 23 or 24, wherein an average distance 
between adjacent transducers increases with increasing distance of said 
transducers from a centre of the array. 

26. The sound system of claim 22, 23, 24 or 25, wherein the one or more digital 
30 signal modifiers is adapted to impose a frequency dependent spatial gain 

window onto the array of transducers. 
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27. The sound system of any one of claims 22 to 26, having digital processors 
adapted to arrange the signal into two or more channels, including the at 
least one rear channel, said channels having different travel lengths to a 
given location, wherein the one or more digital signal modifiers are adapted 

5 to maintain a different beamwidth for each of said two or more channels 

and to impose a frequency dependent spatial gain window onto the array of 
transducers and wherein an average distance between adjacent transducers 
increases with increasing distance of said transducers from a centre of the 
array. 

10 

28. The array, method or sound system of any one of claims 1 to 27, wherein 
the array is a two-dimensional array. 
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FIG. 1 (prior art) 
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FIG . 6B 
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FIG. 7 
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